Proton-proton bremsstrahlung below and above pion-threshold: 
the influence of the A-isobar. 



F. de Jong(i'2'3), K. Nakayama^^'^), T.-S. H. Lee^^) 

Department of Physics and Astronomy, University of Georgia, Athens, GA 30602 

^"^^ Institut fiir Kernphysik, Forschungszentrum Jiilich, 52428 Germany 

, ^ ' Kernfysisch Versneller Instituut, 9747 A A Groningen, the Netherlands 

^ ■ 

Q ! Physics Division, Argonne National Laboratory, Argonne, IL 60439 

Q 

00 



> 



X 



The proton-proton bremsstrahlung is investigated within a coupled-channel 
model with the A degree of freedom. The model is consistent with the A^A^ scattering 



O 

^Sj ' up to 1 GeV and the 7A^A vertex determined in the study of pion photoproduction 

'nI" ' reactions. It is found that the A excitation can significantly improve the agreements 
0\ ; 

^ ' with the pp pp'-/ at Eiab = 280 MeV. Predictions at Eiab = 550 and 800 MeV are 

Q ' presented for future experimental tests. 



I. INTRODUCTION 

The Nucleon-Nucleon bremsstrahlung (NN'-f) reaction has long been considered as a tool 
to investigate the off-shell behaviour of the A^A^ interaction. With the availability of new 
experimental data of pp —>■ ppy reaction |1|,0] , theoretical interest has been revived recently 
0^J^,|,^|[. As a start, it is sufficient to consider only the one-nucleon current to investi- 
gate this reaction since the leading interaction current due to one-pion-exchange vanishes 
identically in the ppy reaction. In a calculation taking into account the relativistic features 
of the one-nucleon current , the sensitivity of the the pp pp~f spin observables to the 
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A^A^ off-shell dynamics was demonstrated. Although the calculated analyzing powers are in 
general in good agreement with the data, the cross-section results are more controversial. In 
Ref. [|I| a normalization factor of 2/3 was proposed to account for the discrepancy with the 
data. However, this discrepancy might also indicate the importance of new mechanisms not 
included in the calculations, thus omitting the need of this rather arbitrary normalization 
factor. In Ref. we have explored qualitatively the mechanisms due to the A degree of 
freedom. The main objective of this work is to present a consistent approach to quantify 
our investigation. This is accomplished by extending the coupled-channel formulation of 
A^A^ scattering developed by Lee and Matsuyama [|l3;|n|] to include the electromagnetic 
interactions introduced by de Jong et al. p. 

The effect of the A excitation on NN'y has been investigated in the past using very crude 
models. Bohannon, Heller and Thompson [|1^ derived the one-pion-exchange A current by 



taking the static nucleon limit of the pion photoproduction amplitude. The A width was 
neglected and hence their model is limited to low energies. They found that the A effect 
is negligibly small in np bremsstrahlung. No calculation of pp bremsstrahlung based on 
their A current has been reported. In contrast to the npy process, where the effect of 
the A-isobar is greatly reduced by cancellations between the different diagrams involved 



due to iso-spin factors, the A-effects may be significant in pp'j reactions. Tiator et al. |T3 
evaluate the contributions of the radiative A-decay in Born-approximation and add these 
incoherently to the nucleonic contributions which were calculated using the Soft Photon 



Approximation (SPA). Although Szyjewicz and Kamal [|14| evaluated all one-pion-exchange 
A contributions and add them coherently, their ppy calculation was carried out using only 
the Born approximation. In our previous work 0, we added the single-scattering A-decay 
diagrams coherently to a state of the art calculation of the nucleonic current contribution 
and found a significant effect. In the present paper we go beyond this in two respects. 
Firstly, we use off-shell T-matrix elements that are generated from a meson-exchange coupled 
A^A^ © A^A © ttNN model which is constrained by the A^A^ scattering up to 1 GeV. The 
quality of the phase-shift fit of this T-matrix is superior to the fit of the T-matrix used in 



Ref. 0. Secondly, we also include the rescattering diagrams induced by the presence of A 
degrees of freedom. We follow the approach of Ref. to account for the relativistic features 
of the and A currents. We will show that, even at energies below the pion threshold, 
the influence of the A on the pp7 reaction is considerable and cannot be neglected in a 
quantitative comparison with the data. 

In addition to giving a more complete description of NN'y reaction at low energies, our 
approach is appropriate for making predictions in the intermediate energy region where 
the pion production through the A excitation becomes crucial. A good understanding of 
the NN'-f reaction at intermediate energies is needed to describe the production of hard- 
photons in intermediate energy heavy- ion collisions p!5|J16[| . In these complicated processes 
the photon provides a clean probe of the reaction dynamics. Furthermore, recent studies of 
dilepton production in proton-nucleus collisions have shown that virtual A^A^ bremsstrahlung 
and A-decay are the dominant reaction mechanisms |lT7| , p!8| . The model presented in this 
paper can be used as the starting point of a microscopic approach to investigate both the 
photon and A productions in relativistic heavy-ion collisions. 

In section II, a coupled-channel formulation for NN'y reaction will be presented. The 
results are presented in section III. Section IV is devoted to discussions of future develop- 
ments. 



II. THE HAMILTONIAN MODEL OF NN AND NN-f INTERACTIONS 

To investigate the effect of the A on the NNj reaction, it is more transparent to employ 
a formulation including the A degree of freedom explicitly. In the first part of this section, 
we will describe such an approach based on an extension of the Hamiltonian model developed 
by Lee and Matsuyama p^ , pJ]] to include the electromagnetic couphngs introduced by de 
Jong et al. 0. We then describe how the procedures introduced in Refs. are used to 
include the "minimum relativity" in the NN'-f calculation. 

The objective of the formulation of Refs. 0,|n|] was to obtain a consistent description 
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of A^A^ scattering from the low energy region {Eiab < 400 MeV) to the intermediate energy 
region where the pion production through the A excitation becomes important. In contrast 



with the other nNN models [|T9|, the formulation of Refs. pO| , pn| was designed to have 
a smooth transition to the usual nonrelativistic potential model of A^A^ interaction. This 
was achieved by using a substraction procedure to define the A^A^ potential in the coupled 
A^A^ © A^A © ttNN space from a chosen A^A^ potential which fits the low energy A^A^ data. 
The resulting model not only maintains the good fit to the low energy A^A^ phase shifts, but 
can also describe the A^A^ data up to about 1 GeV. Starting with such a coupled- channel 
model will make the present study significantly different from all of the previous studies of 
the A effect on the NN'y reaction. In particular, we will be able to make realistic predictions 
in the intermediate energy region where extensive data will soon become available at COSY, 
the proton-cooler ring at the Forschungszentrum Juelich (Germany). 



In the formulation of Refs. |lTO| , pl]| , it is assumed that the Hamiltonian for A^A^ scattering 
can be written in terms of three degrees of freedom: A^, vr, and A. In this work, we 
further assume that the NN'y reaction can be described by adding to this Hamiltonian the 
electromagnetic interaction Vgm introduced in Ref. |^. The model Hamiltonian then takes 
the form 

H = Ho + Hint + Vera, (l) 



where Hq is the sum of relativistic free energy operators, e.g. + for the a = N, 
A and vr degrees of freedom. We will neglect the nonresonant pionic interactions {H2 of 
Eq.(1.2) of Ref. [|ll|]) which are found in Ref. ||lT| to be unimportant in describing A^A^ 



scattering. The A^A direct interaction is also neglected in Refs. [p!0| , pn^ for simplicity. The 

considered hadronic interaction is then of the following form 

2 

Hint = "^[KNAi^) + ^A,7rAf(0] 
1=1 

1 2 

+ o H [VNN,NN{hj) + VNN,NA{hj) + VNA,NN{hj)] (2) 



The derivation of the A^A^ scattering equation from the hadronic Hamiltonian Hq + if, 



int 



can be found in in Ref. ||2^ and was summarized in Ref. |]TT|. For our present purpose, 
we will neglect the much weaker effects due to the A^A scattering (induced by the vertex 
interaction /ittAt^a) and the coupling to the nd channel. The resulting scattering T-matrix 
in the coupled A^A^ © A^A space then takes the following form 

Pnn 

TNN,NNiE) = Vnn,nn{E) + Vnn,nn{E)— — — -^Tnn,nn{E), (3) 

-C/ — riQ + te 

Pnn 

TNA,NNiE) = VnA.Nn[^ + , tt — r~-^^^:^^(-^)] ' ('^) 

rj + Uq + te 
Pnn 

TNN,NAiE) = [1 + Tnn,nn{E)— — — —^]Vnn,na, (5) 

-C/ — -Ho + 

where the effective A^A^ potential is 

yNN,NN{E) = Vnn,NN + ^i^AT, AT AT (6) 

with 

UNN,NNiE) = VnN,NAj^ _ JJ^^^J-^^^J^^ Vna,nn- (7) 

In the above equations, Pnn and Pna are respectively the projection operators for the A^A^ 
and A^A channels. The A self-energy Sa is determined by the hT^N*-*A vertex interaction in 
the presence of a spectator nucleon 

2 p 

^a{E) = V/iA,7rAf(«)7; — . Kna^.^)^ (8) 
i=i h — Hq + le 

where Pt,nn is the projection operator for the ttA^A^ state. The vertex interaction h^^N^-^A 
was determined from fitting the vrA^ P33 phase shifts. The transition potential Vnn<^na 
was taken from the one-pion-exchange model of Niephaus et al. with a monopole form 
factor (A^ — m^)/(g^ + A^) {q denotes the three-momentum transfer). The A^A^ interaction 
is defined by an subtraction of any NN potential which fits the A^A^ phase shifts below the 
pion production threshold. In this work we consider the Paris potential and hence the 
A^A^ interaction in Eq. (|^) is defined by 

VnN,NN = Vpar^s ' U^'^ (E = E,) (9) 



where the substraction energy Es is a parameter. The above definition of the A^A^ interaction 
amounts to removing phenomenologically the two-pion-exchange with an intermediate A'^A 
state from the Paris potential, in order to avoid the double counting of the A^A effect. In 
Refs. [|10],|ni, it was found that by choosing A = 650 (MeV/c) and Es = 10 MeV (the 
laboratory energy of the incident nucleon) the solution of Eq. (0) can best reproduce the 
phase shifts calculated from the Paris potential at energies below about 300 MeV, and also 
can describe the A^A^ phase shifts reasonably well up to about 1 GeV. 

The numerical method for solving the above coupled-channel equations Eqs. ®-® in 
the momentum-space representation was well developed in Refs. [p!OipA| . The calculated 
plane- wave matrix elements of T^j^f^NN^E), T^j^^^a^E) and T/vA,ArAf(-E') are the input to the 
study of the NN'j reaction. The formalism presented above yields T-matrices which are 
a solution of the non-relativistic Schroedinger equation. Thus the expression for the A^A^ 
cross-section calculated with these T-matrices is 



da m 



TNN,NNiE)\^ (10) 



dn An 

where ttt-at denotes the nucleon mass. 

Before we proceed further, it is necessary to define the electromagnetic interaction Vem 
of Eq. (|I]). For the pp pprf reaction, it is sufficient to only consider the one-baryon cur- 
rents since the leading term of the two-body current is absent in the pp^ reaction. Similar 
to the approach developed in Ref. [^] to study the electromagnetic production of pions, 
the one-baryon current operator introduced in Ref. ^ is defined by the matrix elements 
of the Feynman amplitudes calculated from an effective Lagrangian. In this way we in- 
clude the important relativistic spin-correction into our model. The matrix-element of the 
electromagnetic transition potential is then given by: 



{p',kX\Km\P) = J^yrkX\Km\P)J^. (11) 



Where {p' k, X\Vem\p) denotes the Lorentz- invariant matrix element and E^ = \J {p^ + rn?^). 
For the A^ ^ N'j vertex this is [k + p' = p) 



ij>\k\\V,m\p)=u{p')[-ie^^+ mUP) (12) 

where A denotes a photon state with polarization e^. We follow the conventions of Bjorken 



and Drell p4| The nucleon Dirac spinor is denoted as normalized as u{p)u{p) = 1. The 
charge and anomalous magnetic moment of the ith nucleon are denoted respectively by 
and /ij. Similarly, the matrix elements of the 7A and A vertex interactions are 

respectively 

{pA,kMVem\p) = r{pA)T^''Mp) (13) 

and 

(P, kX\Vem\pA) = u{p)K^'r{PA) (14) 

where iP^{Pa) is the Rarita-Schwinger spinor with normalization ip'^{p)'ip^{p) = —1- We 
follow Jones and Scadron [ESl to write the gauge-invariant vertex functions: 



^NA, = Kl + Kl (15) 

with 

Kl = ieG^{^t^- ^,)i^T, (16) 

Kl = ieG2{t,P -k-e- Pk^)r'T,. (17) 

For the decay of a A in a nucleon and photon we have: 

= -K + K- (18) 

In the above expressions fc^ = p^^ — is the photon momentum (defined to be outgoing 
from the vertex) and P = |(p^ + P^)- is third component of the isospin transition matrix 
for coupling an isospin 3/2 to an isospin 1/2 particle. 

The coupling constants Gi and G2 in Eqs. (|T6|) and ( pT[ ) are conventionally determined 
by fitting to the M1+ and El^ multipole data on the photoproduction of pions from nu- 



cleons 1 23, 25. 26, 27, 28]. The values obtained depend on the treatment of the non-resonant 



background contributions. Although this leads to some uncertainty in the values, the pa- 
rameters found in the literature are not too far apart. It has also been shown p7| , |28| that in 
order to accurately reproduce the Ml"*" multipole data on the pion photoproduction around 
the resonance energy, one needs energy-dependent couplings Gi and 6*2 • However, given the 
uncertainty in the coupling constants we ignore this dependence. Bearing in mind that the 
vertex Kj^ gives the dominant contribution, we can classify the various sets of coupling con- 
stants by the magnitude of Gi. The lowest value is found by Nozawa et al. [^]: Gi = 2.024 
(GeV~^) and G2 = —0.851 (GeV~^). Highest values are given by Jones and Scadron 



Gi = 2.68 (GeV-i) and G2 = -1.84 (GeV-^) and by Davidson et al. |6l: d = 2.556 



(GeV-i) and G2 = -1.62 (GeV^^)^ ^ote that the authors of Ref. ||] have a slightly dif 



ferent definition of K2. An alternative way to extract the A^A7 coupling parameters is to 
assume vector- meson dominance. On the NA'y vertex only the isospin-1 vector meson con- 
tributes and the coupling strengths are determined by the ratio of gpNN and gpNA- Using the 
coupling constants of Ref. [p3] this procedure gives Gi = 2.0 (GeV~^) and G2 = (GeV~^), 



comparable with the values from pion photoproduction. 

In the description of A^A^ bremsstrahlung reactions, where the calculations are based on 
A^A^-potential models, one makes use of the Lorentz invariant nature of the various NN'j 
transition amplitudes describing different bremsstrahlung processes. It is then important 
to be able to obtain such invariant amplitudes from potential models that are based on 
non-relativistic approaches, which yield Galilean invariant amplitudes. This is the case of 
the present model given by Eqs. (0)-(|D. In Ref. a procedure is described to construct 
Lorentz-invariant amplitudes from non-relativistic amplitudes. Following this prescription 
we construct an amplitude that transforms covariantly by including the proper E/m factors 
and using the proper relativistic kinematics. This procedure is similar to the introduction 



of "minimal relativity" in earlier studies of relativistic effects on nuclear bound states pO 
However, although the procedure provides us with matrix elements that transform covari- 
antly, our formulation is not fully relativistic. It suffers from the inconsistency commonly 
encountered in constructing realistic meson-exchange A^A^ models. This deficiency is un- 
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avoidable in practice since we do not have a reliable relativistic theory for describing the 
short-range interactions which are clearly beyond the meson-exchange description. We also 
point out that our formalism is equivalent to what one obtains when starting out from a 
field-theoretic point of view and: a) approximates all propagators by their positive energy 
content, and b) applies a three-dimensional Thompson-like reduction to all integrations over 
the four-momentum. 

With the "minimal relativity" prescription, and in first order of the electromagnetic 
coupling, the Lorentz-invariant amplitude of the pp pp'f reaction defined by the model 
Hamiltonian Eq. is 

{pp-f\M{E)\pp) = {pp-f\M'^'\pp) + {pp-f\M'^''\pp) (19) 

where the external, or single-scattering, term is 



(pp7|M^^*bp) = J2 {pin\Vem\a)G^{E){a\fMa,NN{E)\pp) 

a=N,A 

+ {ppITnnME - E^)\a)Go,{E - E^){-fa\Vem\pp) (20) 



a=N,A 



Ey is the photon-energy and we have introduced the and A propagators {E^ = \Jp^ + rn? 



N 

etc.) 

^^^^^ ^ Epp'-E.^-w^'^ip'^py ^^^^ 

where po = W — Ep-^, the energy available for propagation in the presence of the spec- 
tator particle, which has momentum pi. T/vAr,Aro denotes the Lorentz invariant scattering 
amplitude, constructed from TNN,Na by |5| 

{p'\f^N,Na\p)=(^] {p'\TNN,Na\p) (23) 

The rescattering term is 

(pp7|M--|pp) = Y {pp\TnnME - E,)\5l3)G,p{E - E,) 

a,l3,5=N,A 

X {^6\Vem\a)Gaf,{E){6f3\f^fs,NN{E)\pp) 
9 



with 



''E-E,-E^^-E^,,E-E^^-E^,' 

- - "E^' % - % ^ ^P - ^P - iP ^P) (24) 

In the calculation of Eq. ^ we neglect the terms were more than one of the intermediate 
particles is a A. Thus we do not include diagrams with a two-A propagator (which is 
expected to be small due to the mass difference) and the diagrams with a AA7 coupling. 
In the above equations = 1.23 GeV is the physical mass of the A, and ttia is the 



bare mass of the A. Also E^* = yp"^ + and E^ = \Jp'^ + m\. The A propagator we 



use is consistent with the relativistic extension of the A- propagator used in Ref. p9|. The 
specific form of S^(p^,pa) determines the value of the bare A mass due to the requirement 
that the A propagator is resonant at the physical resonance position, = m^. The form of 
Sa is restricted by the P33 phase-shift. However, one finds satisfactory reproductions of the 
P33 phase-shift with rather different forms of S'^. For example, the Bransden-Moorhouse 
parametrization only has an imaginary part of the self-energy, implying m\ = m\. In the 



approach of ter Haar and Malfliet ||29[ one gets = 1.46 GeV. Lee |T0| reports an accurate 
reproduction of the P33 phase-shift over a large energy range for a value = 1.28 GeV. 
As we will show later, the form of S'^ has a definite influence on the magnitude of the A 
decay diagrams. 

The various processes contributing to the NNj reaction defined by Eqs. (pOD-(^4D are 
illustrated in Fig. |l|. They can be classified as: (1) the contributions from the nucleon 
current (Figs. |l|a and |l]b), (2) the contributions from the direct A-decay (Fig. |l|c), and (3) 
the contributions from the A^A rescattering (Fig. |l]d). Note that the A^A^ and A^A T- matrix 
in each diagram is determined by different collision energies since the outgoing photon shares 
the total energy available to the system. This complicates the calculation if we follow the 
conventional approach based on the partial-wave decomposition. Instead, we directly carry 
out the calculation of the diagrams with a A^A7 vertex in a plane-wave basis. The plane- 
wave T-matrices are constructed from the partial- wave solutions of Eqs. (|3)-(ll) by including 
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partial waves up to J = 9. Calculating in the plane- wave basis allows an exact treatment of 



the relativistic features such as that of the vertex interactions defined by Eqs. (|T^)-([T8|). 

With the "minimal relativity" prescription described previously, all matrix-elements of 
the process displayed in Fig. are Lorentz invariant. We therefore are allowed to 
calculate each of them in any convenient frame. The calculation of the contributions from 
the nucleon current (Figs. |l]a,b) is identical to that of Ref. P], except that the A^A^ T- matrix 
is now generated from the coupled-channel equations By setting the coupling term 

f/(i) in Eqs. (g) and (g) to zero, we reproduce the results of Ref. [^]. For the single-scattering 
A-decay diagrams (Fig. |l]c), it is most convenient to do the calculation in the c.m.-frame of 
the initial NN state. Explicitly, from Eq. ^ we have for this diagram 



pi 

X (p',r;',r^|fjvA,iV7v(s = 4E„^')|p, ri, rs). 



with 



(/,r;|f^^^|p',r'/) = u^^Xp")Tf^rrl{p). 
{p\ r'l, r'^\fNN,NA{s)\p, ri, ra) = ip'^^ {p')ur'^{-p')f^^'^^ {s)ur,{p)ur^{-p) 

p^ = 2E^ - E^. (25) 

In these equations r denotes the spin index of the respective spinor, r^^''' is defined in 
Eq. (p!B[), the momenta p,p' and p" are determined by the external kinematics. A similar 
expression is found for the second diagram of Fig. |T]c where the photon is emitted before 
the strong interaction. 

The expressions for the rescattering diagrams with A's are evaluated in the same fashion. 
We see from Eq. (p^ that the calculation involves an integration over the intermediate state. 
Again we use the Lorentz invariance of the ppj amplitude and evaluate the diagram in a 
suitable frame. We find for the rescattering diagram where a nucleon is excited into a A 
and a real photon 
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^2E^ -ko- E^ - Ep_-^ - ^SA(pA^p. _ k')^ 



=E^-ko. (26) 

In this expression k is the photon momentum. The notation is similar to Eq. (^), we 
specified exphcitly all the momenta of the AA^ T-matrix which is not in its cm. frame. Due 
to the presence of the (complex) A self-energy the denominator with has no zero in the 
region of integration. The other denominator does have a zero, evaluating the diagram in 
the cm. frame of the A^A^ T-matrix allows the pole to be treated with a simple subtraction 
method. The other rescattering diagrams with a A intermediate state are calculated in the 
same manner. 

Finally we point out that the rescattering diagram where the A is the spectator has an 
intermediate state which separates in a proton-A"*" and a neutron-A++ state. The photon 
couples to the nucleon, hence this diagram has a contribution depending on the magnetic 
moment of the neutron. This contribution is gauge-invariant, however the contribution of 
the proton-proton-photon vertex is not. Also, since the introduction of A intermediate states 
allows for the exchange of charged mesons in the proton-proton T-matrix, meson-exchange- 
currents (MEC) will be non-zero. This all shows that the inclusion of A intermediate states 
has implications on the gauge-invariance of the model. 



III. RESULTS 

In this section, we will first present results in the low energy region where the experi- 
mental data at Eiab = 280 MeV are available. We then present our predictions at Eiab = 550 
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Mev for the forthcoming experimental test at COSY. We will also consider the kinematic 
region where the A dynamics can be best studied. 

Compared with the previous works in the low energy region {Eiab < about 300 MeV), an 
important feature of our approach is that we calculate the T-matrices Tnn,nn and Tnn^na 
from a coupled-channel model which was obtained by extending the Paris potential to include 
the coupling with the A^A and nNN states. As shown in Refs. |10| , [TT1| and briefly discussed 



in section II, the constructed coupled- channel model is as good as the Paris potential in 
describing the NN data in the low energy region below the pion production threshold. In 
order to illustrate the on- and off-shell differences between teh A^A^ interaction based on the 
Paris potential and our coupled channel approach, we decompose the A^A^ T-matrix in its 
spin-isospin compenents ||3l| 



{p'\Tnn,nn\p) = [(yiPs=o + «2-Ps=i + ia-s^ai + ^2) ■ n 

+ a4^i2(g) + «55i2(Q) + tt65i2(g, Q)]Pt- (27) 



where 



p — p' - p + p' ^ p X p' 
—nQ = 



\p — p'\ \p + p'\ \p X p'\ 

Ps=0 = ^(1 - CTl • ^2), Ps=l = ^(3 + (Ji ■ (72) 

Si2{p) = 3^1 ■pa2-p-ai- (T2, 
3 

Suip,?) = 2 ■ P'^^ ■p + oi-pa2- p') - p' ■ poi ■ 0^2. (28) 

The first two terms in Eq. |2^, proportional to the spin projection operator Ps=o and Ps=i, 
are the central spin-single (S = 0) and spin-triple (S = 1) components, respectively. The 
third term is the spin-orbit component, the fourth and fifth terms are the usual tensor 



components. The last term in Eq. 27 is the off-shell tensor component which vanishes 
identically on-shell as a consequence of time-reversal invariance. Pt in Eq. denotes the 
total isospin projection operator. For proton-proton scattering, it projects out onto the total 
isospin T = 1 subspace. In Fig. □ we compare the A^A^ on-shell interaction. Here, following 
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Ref. we show the angle-averaged magnitude of each spin-isospin compenent with total 
isospin T = 1. The dashed curves are the results for the Paris potential while the solid 
curves stand for the coupled-channel model results. The only notable differences are in the 
Ps=i channel. The small differences seen in Fig. ^are expected since the subtraction method 
defined by Eq. (j^) yields a correction U^^\E) — U^^^Eg) to the Paris potential, which is 
very small at low energies and becomes significant only as the collision energy approaches 
the pion production threshold. 

The main difference between the Paris potential and the coupled-channel model is in 
the off-shell T-matrices which describe the wavefunctions in the interaction region. As an 
example, we compare in Fig. ^the half-off-shell T-matrices calculated from these two models 
at Eiab = 280 MeV as a function of the off-shell momentum p' = \p'\. Significant differences 
can be observed in most of the channels, especially at higher off-shell momenta. However, the 
pp~f reaction at this low energy only probes the low momentum region {p' < 0.36 GeV/c) in 
which the differences between the two half-off-shell T-matrices are much smaller. Moreover, 
the ppy reaction is insensitive to the central Ps=o channel which shows the largest 



difference in Fig. ^. We compare in Fig. ^ the NN'y analyzing powers predicted by the Paris 
potential (dashed curve) and the coupled-channel model (solid curve), for proton angles of 
12.4 and 14.0 degrees at El = 280 MeV. To emphasize the farthest off-shell region, we 
did not include the rescattering mechanisms (Figs. |l|b,d) in this calculation. Obviously, the 
off-shell differences shown in Fig. |^ do not lead to any dramatic effects on the NN'y reaction 
at low energies. We mention that the analyzing power is most sensitive to the tensor and 



central spin-triple components of the A^A^ interaction [p2 |. 

We now turn to analyzing the effect due to the A mechanism which was found in Ref. 
0] to be significant already at E^ = 280 MeV. At certain kinematic conditions, the A 
mechanism can increase the cross-section up to 30 %. But it can decrease, less substantially, 
the cross-section at other kinematic conditions. Clearly there is a very large interference 
effect due to the presence of the A excitation. The A-current, which is almost exclusively 
magnetic, interferes very effectively with the magnetic part of the nuclear current. The latter 
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gives the dominant contribution at this energy. To see how this arises within our model, we 
present in Fig. ^the results at Ei^b = 280 MeV calculated from various combinations of single 
scattering mechanisms (Figs. |l]a,c). We first observe that at this low energy the contribution 
from the A excitation alone (dot-dot-dash curve) is less than 10% of the contribution from 
the nucleon current (dot-dash curve). The A contribution depends only very weakly on the 
photon angle, while the nucleonic contribution shows a large variation. Second, we note 
that the A contribution (Fig. |]c) consists of two different amplitudes: (1) the pre-emission 
amplitude due to the emission of the proton before the strong interaction takes place, (2) 
the post-emission amplitude due to the emission of photon after the A^A^ collision. The 
difference between these two amplitudes is mainly in the 7A^A vertex. The pre-emission 
amplitude is determined by the vertex Y^^^, and the post-emission amplitude by Y^'^^ . 
From the expressions Eqs. (|15|) and ([T8|), we see that the dominant Kjj^ term of the pre- 
emission vertex Y^'^^ has a minus sign relative to that of the post-emission vertex Y^^^. 
Consequently, the pre-emission and post-emission amphtudes tend to have opposite effects 
in interfering with the nucleonic contribution. This is also illustrated in Fig. ^. It is seen that 
when the pre-emission amplitude is added to the nucleonic amplitude, the interference effect 
is destructive at all angles and yields the dashed curve in Fig. |^. On the other hand, the 
interference due to the post-emission of A is constructive. By further adding this amplitude 
to the calculation, the result is shifted from the dashed curve to solid curve. This sensitive 
interference between the nucleonic and A contribution provides an opportunity to test our 
model of the A excitation. 

In calculating the A contributions to the NN'y reaction it is crucial to use a A propagator 
which is consistent with the T-matrices used. As discussed in section II, there are other forms 
of the A self-energy found in the literature. All are constrained by fitting the vrA^ phase-shifts 
in the P33 channel, but results from different formulations of vrA^ scattering. To show how the 
form of the A self-energy (in casu the A propagator) affects the results, we show in Fig. ^ a 
results from our calculation (dashed curve) using Eq. (^ to determine the A self-energy and 



a results (dotted curve) using the self-energy of ter Haar and Malfliet The contributions 
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of the nucleonic current are represented by the sohd curve. In all calculations the same T- 
matrices based on the Paris potential including A intermediate states were used. Although 
both propagators fit the P33 phase-shift, the rather large difference between the results 
show the need to use a consistent A propagator. Consequently, all of the results presented 
hereafter are obtained using the self-energy of Eq. in calculating the A propagator Eq. 

The inclusion of the A rescattering diagrams (Fig. |l]d) is one of the main new features 
of this work. It provides an additional contribution due to the A excitation. In Fig. 0, we 
illustrate its effects at various laboratory energies in changing the calculated differential cross 
sections and analyzing powers. The contribution from the nucleonic currents alone is denoted 
by the dotted curves. By adding the A single scattering mechanism (Fig. |I]c), we obtain 
the dashed curves. The solid curves are obtained when the A rescattering mechanisms (Fig. 
|l]d) are also included. At lower energies the effect of the A-rescattering on the analyzing 
power is comparable to that of the A single-scattering mechanism. The contribution of 
the A-rescattering to the cross-section is about one half of the single-scattering mechanism. 
At higher energies the rescattering diagrams are less significant. We have found that the 
largest A effect is due to the rescattering diagram where the A is a spectator (the rightmost 
diagram in Fig. [l|d). The calculation of this mechanism involves an integration over two 
half-off-shell Tj^n^na matrices and hence is most sensitive to the short-range behaviour 
of the transition potential V^j^f^jy^ potential used in solving the coupled equations Eqs. 
(3)"(i)- the construction of Refs. the model of Vatat^-^ata derived by Niephaus et 

al. ||21| is used. Because of the use of a different regularization of the short-range part of 



the one-pion exchange, this model is rather different from the model employed in the NN'y 
calculation of Ref. 0. Consequently, the resulting effect of A rescattering is rather different. 
The A rescattering effect calculated using the Tattv^tva matrices employed in Ref. is much 
smaller in changing the differential cross sections. This is also the reason why the structure 
of the A-contributions to the analyzing power at 280 MeV in the present work differs from 
the results found in Ref. f^. We hope that the new high- precision measurements scheduled 
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at the CELSIUS facility in Uppsala (Sweden) and the COSY ring in Juelich (Germany) will 
shed more light on this issue. 

In Figs. P and ^, we compare our results with the experimental data at Eiab = 280 MeV. 
All mechanisms in Fig. |1| are included in the calculations. The scattering T-matrices and 
the A self energy needed in evaluating the ■jNN amplitudes defined by Eqs. ([TP|)-(P^ are 



generated from the model of Refs. |[T^JTT|, as discussed in section II. Since there are still 



some uncertainties in determining from the pion photoproduction reaction the values of Gi 
and G2 of the 7A^A vertex, we performed calculations for two sets of coupling constants. 
The dashed curves are calculated from using the highest value of Gi = 2.68 (GeV^^) and 
G2 = —1.84 (GeV~^) as determined by Jones and Scadron ||25|. As for their application in 



ppy calculations these values are very close to the ones found in a recent analysis of Lee, 
who finds slightly higher values H: Gi = 2.89 (GeV-^) and G2 = -2.18 (GeV^^)^ jhe 



solid curves are from using the lowest value of Gi = 2.0 (GeV~^), G2 = 0.0 (GeV~^) as 
predicted by the vector dominance model. As a reference we also include the results (dotted 
curves) calculated from only the nucleonic contributions (Figs. ^,b). As can be seen in Fig. 
^, the A excitation mechanisms increase the cross-section in most of the range of kinematics 
considered. The A-decay diagrams (Fig. |I|c) are responsible for the largest part of the 
increase. However, the nucleonic contributions are still dominant at this low energy and the 
differences due to the use of two different sets of 7A^A coupling constants are not particularly 
large. Nevertheless, the A effects clearly significantly bring the theoretical values closer to 
the data of the differential cross sections. 

The A effects also improve the agreements with the data of the analyzing powers, as 
seen in Fig. ^. The agreement with the data is very good in most of the kinematic regions 
considered. Again, the differences between the results using two different sets of Gi and G2 
are rather small. Clearly, it is necessary to consider the higher energy regions in order to 
have a critical test of our model of the A excitation and pin down the values of Gi and G2 
of the 7A^A vertex. 

To illustrate the increasing importance of the A mechanisms as energy increases, we show 
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in Fig. the inclusive photon production cross-section as a function of photon energy at 
photon angle 9 = 90° and Eiab = 280, 550 and 800 MeV. At Eiab = 280 MeV we find a small 
decrease in the cross-section for photon energies smaller than 80 MeV and a minor increase at 
higher photon energies. Although the cross sections are clearly dominated by the nucleonic 
contribution, the A effect is still significant in influencing the exclusive cross sections as 
displayed in Figs. | and p. At higher laboratory energies we clearly see the fc-dependence of 
the A-decay contributions, resulting in peaks at the high end of the kinematically allowed 
photon energies. The cross-section near the A peak is doubled by the A mechanisms at 
Eiab = 550 MeV, and by a factor of about 5 at Eiab = 800 MeV. The large enhancement of 
the A peak is mainly due to the resonant behavior of the A propagator Eq. ( p2D and the 
dependence on the photon momentum of the vertex Eqs. (|15|,|18D. This is in contrast to the 



results at Eiab = 280 where the effect is due to strong interference of a small A-contribution 
with the much larger nucleonic current contribution. 

In the figure for Eiab = 800 MeV we also plot the cumulative contributions of the partial 
waves for J < 1 to J < 4. The largest contributions are from the J = 3 partial wave(^-F3). 
It is interesting to point out that the same large contribution from the J = 3 partial wave 
is also found in the study of pp — > pnir'^ reaction (Table I and III of Ref. To have a 



unified description of both the pionic and electromagnetic excitation in pp scattering, it will 
be interesting to also have the NN'j data at 800 MeV. 

For future experimental tests of our model of the A excitation, we present our predictions 
of the exclusive cross sections and analyzing powers at Eiab = 550 MeV (Fig. |ll]) and 800 
MeV (Fig. |1^). The calculations are identical to that of Figs. | and ^ The dotted curves are 
from calculations including only the nucleonic contributions. The solid and dashed curves are 
respectively calculated by using the 7A^A coupling constants of vector dominance and Jones 
and Scadron |^ . As seen in both Figs. |Tl| and including the A contributions drastically 
changes the differential cross sections. Hopefully, the forthcoming COSY experiment will 
have enough accuracy to distinguish the solid and dotted curves. This will help narrow 
down the values of Gi and G2 of the 7A^A vertex. Rather surprisingly, the inclusion of the 



A-decay diagrams has much less effects on the calculated analyzing powers. This suggests 
that the A-decay diagrams have at these energies and kinematics a similar spin-structure 
(to the extent this is measured by the analyzing power) as the nucleonic contributions. 

Finally, in regard to ppj experiments at high energies we should mention the efforts 
of Ref. These authors performed a pp'y experiment at Eiab = 730MeV. However, 



due to the experimental set-up they probe a kinematical region which has relatively low 
photon-energies. In the coplanar geometries for which we presented our calculations the 
photon-energy in the initial A^A^ cm. frame ranges from 200 - 350 MeV. In the experiment 
of Ref. the photon-energy is below 150 MeV. The data from Ref. show a rise at the 



high-energy end of the kinematically reachable photon-energies, possibly indicating an onset 
of the A effect. However, as can be seen from Fig. |1^, the A effect reaches its maximum 
at much larger photon energies. Together with the large errors this makes the experiment 
unsuitable for discriminating between the various sets of A^A7 coupling constants. 

IV. SUMMARY AND DISCUSSION 

In this work we have developed an approach to investigate the pp bremsstralung at both 
low and intermediate energies. It is based on an extension of the coupled A^A^ © A^A © 



TcNN model of Refs. [|13,|TT|] to include the electromagnetic coupling with the A^ and A 
currents. The hadronic part of the model is consistent with A^A^ scattering up to 1 GeV. 
The 7A^A coupling is determined in the study of pion photo-production in the P33 channel. 
The relativistic features of the one-baryon current matrix elements are treated exactly by 
performing the calculation directly in momentum space. In addition, the A rescattering 
contribution (Fig. |I|d) is evaluated first time in this field. 

We have shown that the A contribution can interfere strongly with the nucleonic contri- 
bution. As seen in Figs. ^ and ^ the A contribution can significantly improve the agreement 
with the data even at low energies below pion production threshold. Given the uncertainty 
in the normalization of the data, it remains to be seen if the remaining discrepancies with 
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the data are due to these normahzation problems or that higher order interaction currents 
such as those due to the p ^ vr 7 and ^ tt 7 couphngs can account for these discrepancies. 
A quahtative study of these effects in conjuction with the A effect can be found in Ref. |^ . 

At intermediate energies, the A excitation dominates the pp —>■ pp'-f cross sections, as 
displayed in Figs. 0-[T^. Our predictions should be reasonable in guiding the experimental 
efforts. However, the employed coupled- channel model was constructed by fitting the NN 
phase shift data back in 1981 This has to be improved since the phase shift data have 
been changed significantly recently |3^. Another important step is to extend the present 
approach to study np bremsstralung at intermediate energies. This requires the inclusion 
of exchange currents within the coupled-channel model. Our effort in this direction will be 
published elsewhere. 
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FIGURES 

FIG. 1. Diagrams included in the calculation (a single line denotes a nucleon, a double line 
a A intermediate state): single-scattering diagrams with Tnn-nn (a), rescattering diagrams 
with Tnn-nn (b), single-scattering diagrams with Tnn-na (c) and rescattering diagrams with 
Tnn-na (d). 

FIG. 2. Magnitude of the spin-isospin components of the on-shell interaction, averaged over 
the scattering angle as a function of the laboratory energy. The solid line stands for the T-matrix 
including A-intermediate states, the dashed line represents the Paris T-matrix. 

FIG. 3. Magnitude of the spin-isospin components of the half off-shell interaction, averaged over 
the scattering angle as a function of the off-shell momentum at a beam energy Eiab = 280 MeV. 
The solid line stands for the T-matrix including A-intermediate states, the dashed line represents 
the Paris T-matrix 

FIG. 4. Coplanar analyzing power for proton-scattering angles 6i = 12.4 and 02 = 14.0 as a 
function of the photon angle, Eiab = 280 MeV. The full line is the result (including the nucleon 
single scattering diagram. Fig. ||a) calculated with the T-matrix with A intermediate states, the 
dashed line is obtained using the Paris T-matrix. 

FIG. 5. The various contributions to the total cross-section, in a coplanar geometry with pro- 
ton-emission angles 9i = 27.8 and 82 = 28.0, Eiab = 280 MeV. No rescattering contributions were 
included. The dash-dotted line is the result with only the nucleon contributions (Fig. [l|a), the full 
line stands for the full model (Fig. |l|a,c) taking the high NA'y coupling set {Gi = 2.68 (GeV~^), 
G2 = —1.84 (GeV^^)). The dashed line stands for the nucleonic contributions plus the pre-emission 
diagrams. The dot-dot-dash line is the result with only A-decay diagrams (Fig. ^). 
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FIG. 6. The effect of the choice of the A propagator, in a coplanar geometry with pro- 
ton-emission angles 9i = 27.8 and 62 = 28.0, Eiab = 280 MeV. The dotted hne is the result 
with only the nucleon contributions (Fig. ||a), the full line is the result, including all diagrams, cal- 
culated with the A self-energy of ter Haar and Malfliet , the dashed line is the same calculation 
but using the A self-energy of Eq. ^ In these calculations we took Gi = 2.0 (GeV^^) and G2 = 
(GeV-2). 

FIG. 7. The effect of the A rescattering diagrams on both the cross-section and analyzing power 
in a coplanar geometry, at various laboratory energies. The full line is the result including all dia- 
grams, the dashed line is calculated with the nucleonic diagrams plus the single-scattering A-decay 
diagrams (Figs. |l]a,b,c) and the dotted line is calculated with only the nucleonic contributions (Fig. 
|a,b). In these calculations we took d = 2.0 (GeV"^) and G2 = (GeV'^). 

FIG. 8. Coplanar geometry pp^ cross-section in the laboratory frame as a function of the 
photon-emission angle at an incident energy of Eiah = 280 MeV for various proton scattering 
angles 9i and 62 ■ All diagrams were calculated using the NN and A T-matrices of Lee . The 
dotted line corresponds to the calculation which only incorporates the nucleonic diagrams (Figs. 
||a,b). The results including A-decay diagrams are represented by the solid and the dashed lines. 
These correspond to the two choices of the A^A7 coupling constants: the full line stands for the 
results calculated with d = 2.0 (GeV^^) and G2 = (GeV^^) ^nd the dashed one for d = 2.68 
(GeV~^) and G2 = —1.84 (GeV~^). The data do not contain the arbitrary normalization factor 
of 2/3. 

FIG. 9. Same as Fig. |8| for the analyzing power. The data are from Ref. [Q] and have been 
multiplied by a factor of —1. 



25 



FIG. 10. Inclusive cross-section at photon-angle 90° at various laboratory energies as a function 
of photon-momentum. The dashed line is the result with only the nucleonic contributions (Fig. 
||a,b), the sohd hue is the full result including all diagrams calculated with d = 2.0 (GeV^^) and 
G2 = (GeV~^). In the plot for Eiab = 800 we also show the cumulative contributions of the 
partial waves: J < 1 (dotted), J < 2 (dash-dot), J < 3 (dash-dot-dot) and J < 4 (short dash). 
The full result is calculated with J < 9. 



FIG. 11. Same as Figs. ^ and ^ but at an incident energy of Eiab = 550 MeV. 



FIG. 12. Same as Figs. ^ and ^ but at an incident energy of Eiab = 800 MeV. 
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